NASA TM X-183 


S3 13 8 97 

NASA TM X-183 


TECHNICAL MEMORANDUM 

X - 183 

USE OF EXPERIMENTAL STEADY-FLOW AERODYNAMIC PARAMETERS 

IN THE CALCULATION OF FLUTTER CHARACTERISTICS FOR 

FINITE -SPAN SWEPT OR UNSWEPT WINGS AT 

SUBSONIC , TRANSONIC , AND 

SUPERSONIC SPEEDS 

By E. Carson Yates, Jr. 

Langley Research Center 
Langley Field, Va. 

Declassified February 8 , 1963 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON January I960 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL MEMORANDUM X-I 83 

USE OF EXPERIMENTAL STEADY -FLOW AERODYNAMIC PARAMETERS 
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SUMMARY 


Flutter calculations for several swept and unswept wings through 
the transonic speed range have been made by the modified strip-theory 
method of NACA RM L 57 LIO, which employs steady-flow aerodynamic param- 
eters for the undeformed wing. Experimentally determined distributions 
of steady-flow aerodynamic parameters for the undeformed wings were used 
in the present calculations. Comparisons of these calculated results 
with experimental flutter data and with calculations previously made by 
using linearized-theory aerodynamic parameters indicate that the method 
employed gives accurate flutter results for swept wings at subsonic, 
transonic, and supersonic speeds. However, since this method of flutter 
calculation is not applicable when the Mach number component normal to 
the leading edge is near 1.0, it appears that the transonic flutter 
characteristics of unswept wings cannot be calculated by the method as 
given in NACA RM L 57 LIO. An attempt was made to remove this limitation 
by using experimental values of two-dimensional lift-curve slope (instead 
of theoretical values) in the calculation of circulation functions when 
the Mach number component normal to the leading edge was near 1.0. 
Applying this procedure to two unswept wings removed the spurious 
asymptotic rise of flutter speed near Mach number 1.0, but it did not 
result in very close agreement between calculated and experimental flut- 
ter speeds . 

The possibility of obtaining estimates of flutter characteristics 
for unswept wings by use of total aerodynamic parameters has also been 
examined. Although no supersonic calculations were made using total 
aerodynamic parameters, it appears that the procedure outlined gives 
reasonable estimates of flutter characteristics at subsonic speeds. 
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INTRODUCTION 


Reference 1 presented a strip-theory type of flutter calculation pro- 
cedure for finite-span swept and unswept wings based on spanwise distribu- 
tions of lift and pitching moment derived from distributions of aerody- 
namic parameters associated with the undeformed wing in steady flow. 1 
Subsonic and supersonic flutter characteristics for several wings were 
calculated by this modal analysis method and compared with experimental 
flutter data. However, all the calculated results shown in reference 1 
were obtained by the use of steady-flow aerodynamic parameters calculated 
from linearized theory, and hence no calculations were shown for Mach 
numbers near 1.0. Since minimum flutter speeds usually occur for Mach 
numbers near 1.0, it was not possible to evaluate the minimum flutter 
speed by using these theoretical steady-flow aerodynamic parameters. 

Some linearized-theory methods exist for calculating three-dimensional 
oscillating aerodynamic loads at sonic speed (for example, ref. 2). How- 
ever, the application of such calculations to finite-thickness wings which 
involve mixed-flow regions is open to question. The reliability of these 
methods has not been proved in actual use . 

The present report presents the results of flutter calculation for 
several wings through the transonic range, made by the method of refer- 
ence 1, but employing experimentally determined distributions of steady- 
flow aerodynamic parameters obtained from wind-tunnel and flight tests. 

It may be noted that use of experimental steady-flow aerodynamic param- 
eters introduces into the flutter calculation some nonlinear aerodynamic 
effects such as those of finite thickness and viscosity. The magnitudes 
of these steady-flow nonlinear effects are believed to be approximately 
correct for the oscillating wing as long as the frequency is not too high. 

The possibility of obtaining estimates of flutter characteristics for 
unswept wings by use of experimental total aerodynamic parameters is also 
examined herein. 

In the discussion of the limitations of the flutter calculation pro- 
cedure in reference 1, it was pointed out that because a two-dimensional 
lift-curve slope appeared in the expressions for the circulation functions 
used, the method as presented was not applicable when the Mach number 

^In the method of reference 1 spanwise distributions of steady-flow 
section lift-curve slope and local aerodynamic center for the undeformed 
wing are used in conjunction with the "effective 1 * angle -of -attack distri- 
bution resulting from the assumed vibration modes in order to obtain 
values of section lift and pitching moment. Circulation functions modi- 
fied on the basis of loadings for two-dimensional airfoils oscillating in 
compressible flow are employed to account for the effects of oscillatory 
motion on the magnitudes and phase angles of the lift and moment vectors. 
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component normal to the wing leading edge was very close to 1.0. (The 
theoretical two-dimensional lift-curve slope, employed in the calcula- 
tion of circulation functions in reference 1, approaches infinity as the 
Mach number component normal to the wing leading edge approaches 1.0.) 
The possibility of removing this limitation by using experimental values 
of two-dimensional lift-curve slope in the circulation functions is 
examined in the present report. 

Flutter characteristics calculated by using experimentally deter- 
mined distributions of static aerodynamic parameters are herein com- 
pared with experimental flutter data, and with the flutter character- 
istics calculated for the same wings in reference 1 from the theoretical 
distributions of static aerodynamic parameters. Uncoupled vibration 
modes are employed in all flutter calculations. Calculations are shown 
for three swept wings of aspect ratio 4.0, taper ratio 0 . 6 , quarter- 
chord sweep angle 45°, and local center -of -gravity positions at approxi- 
mately 34 percent chord, 46 percent chord, and 58 percent chord. Calcu- 
lations are shown for two unswept wings of aspect ratio 4 . 0 , taper ratio 
0 . 6 , quarter-chord sweep angle 0 , and local center -of -gravity positions 
at approximately 45 percent chord and 59 percent chord. Finally, some 
flutter characteristics calculated by using experimentally determined 
total static aerodynamic parameters are shown for the two unswept wings 
mentioned previously and for a wing of aspect ratio 4.0, taper ratio 1.0, 
sweep angle 0, and local center of gravity at 50 percent chord. 


SYMBOLS 


A 


aspect ratio of full wing including fuselage intercept 


a 


nondimens ional distance from midchord to elastic axis measured 
perpendicular to elastic axis, positive rearward, fraction of 
semichord b 


ac 


nondimens ional distance from leading edge to local aerodynamic 
center (for steady flow) measured streamwise, fraction of 
streamwise chord 


ac 


■n 


nondimens ional distance from midchord to local aerodynamic 
center (for steady flow) measured perpendicular to elastic 
axis, positive rearward, fraction of semichord b, see 
equation (2) 


ac N 


value obtained by applying equation (2) to total ac value for 
the wing. Total ac is defined as the nondimens ional dis- 
tance from leading edge of mean aerodynamic chord to wing 
aerodynamic-center position (for steady flow) measured 
streamwise, fraction of mean aerodynamic chord. 


4 


b semichord of wing measured perpendicular to elastic axis 

b r semichord of wing measured perpendicular to elastic axis at 

spanwise reference station t) = 0.75 

C complex circulation function, F + iG 

local lift-curve slope for a streamwise section in steady flow 

Cj local lift-curve slope for a section perpendicular to elastic 

a ' n axis in steady flow 

total lift-curve slope for wing in steady flow 
a 


F real part of complex circulation function C 

G imaginary part of complex circulation function C 

k^ reduced frequency based on spanwise reference station (t] = O.75) 

b-pO) 

and on velocity component normal to elastic axis, =■ 

V cos Aga 

M stream Mach number 

Mle Mach number component normal to the leading edge 

V flutter speed, measured parallel to free stream (experimental 

values or values calculated by the method of ref. l) 

V-p calculated reference flutter speed obtained by using C7 = 2 rt. 

n ^a.n 

ac n = - and C = Fj + iGj 

P LE Jl M LE 2 - A 

A sweep angle, positive for sweepback 

A taper ratio of full wing including fuselage intercept 

Tj nondimensional coordinate (either spanwise or along elastic 

axis) measured from wing root, fraction of exposed panel 
span or fraction of wing length 

p air density 


5 
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a) circular frequency of vibration 

% circular frequency of first uncoupled torsional vibration 

mode of wing measured about elastic axis 

Subscripts : 

c/4 quantities associated with wing quarter chord 

c/2 quantities associated with wing midchord 

•4c quantities associated with wing 40 percent chord 

ea quantities associated with wing elastic axis 

C circulation functions obtained from oscillatory aerodynamic 

coefficients for a two-dimensional wing in compressible flow 

I circulation functions for two-dimensional incompressible flow 

PROCEDURE FOR USING EXPERIMENTAL STEADY-FLOW AERODYNAMIC 
PARAMETERS IN THE FLUTTER CALCULATIONS 


Since the experimental steady-flow aerodynamic parameters used 
herein were obtained from several sources, and since these parameters 
are applied to several wings, table I has been prepared as a guide to 
the flutter calculations presented in this report. 


Spanwise Distributions of Lift-Curve Slope 
and Aerodynamic Center 

Wings with 4^° sweepback .- Experimental spanwise distributions of 

steady-flow lift-curve slope C^ and aerodynamic center ac for the 

a 

plan form of A = 4.0, A = 0.6, A c ^ = 45° were obtained from the 

data of references 3 and 4 over the Mach number range from 0.6 to 1.2. 
It should be noted that the wing of references 3 and 4 was 6 percent 
thick (NACA 65 AOO 6 airfoil) while the experimental flutter data of 
references 5 and 6 with which calculated flutter characteristics are to 
be compared, were obtained with wings that were 4 percent thick (NACA 
65 AOO 4 airfoil) . No correction has been applied to the steady-flow 
aerodynamic parameters to account for this difference of thickness. 
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Since the modified strip analysis of reference 1 utilizes aerody- 
namic and structural quantities associated with wing strips normal to 
the wing elastic axis, equations (B6) and (B 7 ) of reference 1 were used 
to convert and ac (associated with the streamwise direction) to 

CL 

C 7 and. ac n (associated with the direction normal to the elastic 
a,n 

axis). (See solid curves of figs. 1 and 2.) For convenience, these 
geometrical relations are repeated here in the notation of the present 
report : 



Crossplots of these 
her M for several 


values of C 7 and ac n against stream Mach num- 
a,n 11 

spanwise stations are shown in figure 3 - 


Direct application of equations (l) and (2) for supersonic speeds 
immediately raises a question with regard to the q station at which 
the tip begins to affect and ac n . As shown in the accompanying 

sketch, 



Mach line 
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a consideration only of strips normal to the elastic axis, indicates 
that the tip influences strips only outboard of point S. However, the 

experimental distributions of Cj and ac obtained from the data of 

oc 

references 3 and 4 are associated with the streamwise direction, and 
hence the tip affects all strips outboard of point N. In order to 
adapt the streamwise data near the tip to the strip orientation employed 
in the flutter calculation, the Cj and ac n distributions obtained 

LL y II 

from equations (l) and (2) were cut at station N, and the portions out- 
board of that point were affinely compressed into the q range outboard 
of point S by the relation 

1 - Wrectea = (% S 1 S l) (3) 

The and ac n values inboard of point N were, of course, left 

unchanged. In the gap thus created between points N and S the C7 

t a,n 

and ac n curves were closed by a faired line which as nearly as possible 

represented an extension of the inboard values. Thus tip effects are 

confined to wing sections outboard of point S, and the faired loading 

between sections N and S logically represents a continuation of the 

inboard loading and is unaffected by the tip. Values of C 7 and ac,. 

‘■a,n 11 

obtained in this manner are given by the dashed curves of figs. 1 and 2. 

Figures 1 and 2 indicate that compared to values of C 7 and ac„ 

6 a,n n 

obtained directly from equations (l) and (2), the result of the above 
alteration is to increase the lift-curve slope and to move the aerody- 
namic center rearward in the tip region. These changes will in general 
oppose each other in their effect on the calculated flutter speed. For 
supersonic Mach numbers flutter calculations have been made both with the 
Cl <x n and aCn values obtained directly from equations (l) and (2) 

(solid curves of figs. 1 and 2) and with C 7 and ac^ values altered 

t a,n n 

as described above (dashed curves of figs. 1 and 2). 

From the preceding sketch it also appears that a strict integration 
of lifting pressure in the direction normal to the elastic axis will 
yield finite loads as far outboard as point T. However, in accordance 
with the strip-analysis method of reference 1, no loads are considered 
outboard of q = 1.0. 

The wing of references 3 and 4 was statically tested in the presence 
of a fuselage with maximum radius equal to 13.9 percent of the wing semi- 
span. The flutter data with which the calculated flutter characteristics 
for the 45° wings are to be compared (refs. 5 and 6) were obtained in the 
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presence of a fuselage with radius equal to 21.9 percent of the wing 

semispan. Therefore, only the values of C 7 and ac n outboard of 

a,n Li 

the 21.9-percent-semispan station were used in the flutter calculations. 
That is, t) = 0 corresponds to 21.9-percent semispan. The effects of 

the different fuselage sizes on the distributions of C 7 and ac n 

^a,n 

are generally confined to the neighborhood of the wing root. Such 
effects should therefore have negligible influence on the calculated 
flutter characteristics . 

Unswept wings with taper ratio of O.6 .- Experimental spanwise dis- 
tributions of steady-flow lift-curve slope Cj and aerodynamic center 

ac for the plan form of A = 4.0, Tv = 0.6, A c /4 = 0 were obtained 

from the data of reference 7 over the Mach number range from 0.6 to 
1.05 and from unpublished flight test data for the X-1E airplane at 
M = 1.4l. The wing of reference 7 had A =4.0, A = O. 5 , ^ c /2 = 

and hence had a plan form which was slightly different from that of the 
flutter-tested wings (ref. 5)- However, the wing of reference 7 a s well 
as the flutter-tested wings had NACA 65 AOO 4 airfoil sections. The wing 
of the X-1E airplane had A = 4.0, A = 0.5, A^ c =0, and a modified 

NACA 64A004 airfoil, and thus also differed slightly in plan form and 
airfoil from the flutter-tested wings. No correction was applied to the 
steady-flow aerodynamic parameters to account for these small differ- 
ences, and the distributions of C 7 and ac^ obtained from the data 

of reference 7 and from the unpublished data for the X-1E airplane were 
used directly in the flutter calculations . 

For the unswept wings the effect of the wing tip on the steady-flow 
aerodynamic parameters for each wing strip is properly accounted for 
without application of a tip correction of the type previously discussed 
for swept wings. For unswept wings the freestream direction and the 
direction normal to the elastic axis essentially coincide so that 


= Ci 


a.n 


a 


and 


ac n = 2ac - 1 


The wing of reference 7 was statically tested in the presence of a 
fuselage with maximum radius equal to 12.35 percent of the wing semispan, 
and the maximum fuselage radius of the X-1E airplane was 20.5 percent of 


\ j \ o \ •£- 1 
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the wing semispan. The flutter-test wings of reference 5 "were tested in 
the presence of a fuselage with radius equal to 21.9 percent of wing 
semispan. The treatment of these different fuselage sizes is the same 
as that described previously for the 45° swept wings . 

The spanwise distributions of C7 and ac n obtained from the 

6 a,n n 

data of reference 7 are shown in figure 4. Crossplots of these values 

C7 and ac n against stream Mach number for three spanwise s ta- 
co, n 

tions are shown in figure 5. Spanwise distributions of C 7 and ac n 
J ^a,n n 

t obtained from flight-test data for the X-1E airplane at M = 1.4l are 

) shown in figure 6 and are compared with values from reference 1 calcu- 

j la ted from linearized supersonic -flow theory. 


Total Lift-Curve Slope and Aerodynamic 
Center for Unswept Wings 

The feasibility of using experimental total lift-curve slope C T 

MX 

and aerodynamic center ac^ in flutter calculations has been investi- 
gated for three unswept wings. It has been observed for a number of 
unswept wings that the spanwise distribution of C^ J Cj^ does not vary 

greatly as stream Mach number M increases from subsonic to transonic 
to low supersonic. Furthermore, the spanwise distributions of experi- 
mental have been found to be close to those calculated by the 

method of reference 8. Therefore, the total lift-curve slope Ct is 

a 

introduced into the flutter calculation as a multiplying factor applied 
to the calculated spanwise distribution of C ^ jC ^ obtained from the 

method of reference 8 for M = O.75. Thus 



In the absence of a reliable method for estimating spanwise distribution 
of aerodynamic center at transonic speeds it is assumed in all flutter 
calculations involving total aerodynamic parameters that ac is constant 
across the span and that 


= ac 


ac 


total 


(5) 
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where ac total referred, to the mean aerodynamic chord. Although 

equation ( 5 ) may be reasonably valid for unswept wings at subsonic 
speeds, it is realized that equation ( 5 ) is not a good approximation 
to the loading characteristics of a wing at transonic and supersonic 
speeds. It does, however, permit the introduction of the qualitative 
variation of aerodynamic center with Mach number. The appropriateness 
of equation ( 5 ) will, of course, be determined by the results obtained 
by using it. 

Unswept wings with taper ratio of 0.6 .- For the plan form of 
A = 4 .0, A = 0.6, A = 0, experimental values of and ac^ in 

the Mach number range from 0.6 to 1.15 were obtained from the data of 
reference 9* These values are shown in figure 7 together with the span- 
wise distribution of Cj for M = 0.75 calculated by the method of 


reference 8. 


Unswept wings with taper ratio of 1.0 .- For the wing of A = 4.0, 

A = 1.0, = 0, experimental values of and acjj were obtained 

from reference 10 in the Mach number range from 0.6 to 1.10 and from ref- 
erence 11 in the Mach number range from 0.4 to I.05. The wing of refer- 
ence 10 had an NACA 65AOO4 airfoil, and that of reference 11 had an NACA 
63A004 airfoil. The wing for which the flutter-test data were given in 
reference 12 had a 4-percent-thick hexagonal airfoil. The and acjy 

values obtained from the data of references 10 and 11 are shown in fig- 
ure 8 together with the spanwise distribution of Cj for M = 0.75 cal- 

CL 

culated by the method of reference 8. 


Use of Experimental Two-Dimensional Lift-Curve Slope in the 


In the flutter calculation procedure of reference 1 the expression 


a 


Calculation of Circulation Functions 



was used for the complex circulation function. In this expression 


VJI ON t -1 
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where l^, Z^, l £ are aerodynamic coefficients given in refer- 

ence 13 for two-dimensional airfoils oscillating in compressible flow, 

and C7 is the two-dimensional steady-flow lift-curve slope associ- 

a,n 

ated with the Mach number component M^g normal to the leading -edge . 


In the calculations of reference 1 linearized-theory values were 

used for the C 7 in equation (6). That is, 
l a,n * 


= — for Mr* < 1 
l a,n P^ ^ 


= J_ 
Ple 


for Mlj; > 1 


(7) 


where p^ = |/| M LE 2 " ^ • If equation ( 7 ) is used in equation (6), it 
can be seen that Fq = 0 for M T ^ = 1. Consequently, a range of Mach 
number in the vicinity of M^ = 1 is inaccessible to the method as 

presented in reference 1. The possibility of removing this limitation 

by using experimental values of two-dimensional steady-flow C 7 in 

6 a,n 

equation (6) is investigated herein for the two unswept wings with taper 
ratio of 0.6. 


Experimental values of two-dimensional C 2 for the NACA 6^A00k 


airfoil were obtained from reference 14 over the Mach number range from 
0.8 to I. 25 . These values are reproduced here in figure 9* Figure 9 

C l r 

also shows the variation with Mach number of two-dimensional 


‘'a, exp 




Values of F r 


^a,theor 

calculated by use of the theoretical values of equa- 


tion ( 7 ) can be converted to the experimental 


Cj base (for the NACA 
00 


6^A00k airfoil) by dividing the Fq by these values of 


Cl 


a-, exp 


Cl 


a, the or 
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RESULTS AND DISCUSSION 


Presentation of Results 

The wing configurations for which flutter calculations were made 
are summarized in table I. Table I also indicates the sources of the 
experimental steady-flow aerodynamic parameters used. Flutter charac- 
teristics have been calculated by use of experimentally determined 
spanwiSe distributions of steady-flow lift-curve slope and aerodynamic 
center (figs. 1 to 6) for five wings of aspect ratio 4.0, taper ratio 0.6, 
and quarter-chord sweep angles 0° and 45° . Flutter characteristics based 
on experimental total aerodynamic parameters (figs. 7 and 8) have been 
calculated for three unswept wings of aspect ratio 4.0 and taper ratio 
0.6 and 1.0. Exploratory flutter calculations employing experimental two- 
dimensional steady-flow lift-curve slopes in the circulation functions 
for M t ^ near 1 have been made for the two unswept wings of taper ratio 

0.6. Structural data as well as experimental flutter data for all of 
these wings were obtained from references 1, 5, 6, and 12. 

Since calculated flutter characteristics appear not to be very sensi- 
tive to slight changes in mode shapes, and since the wings of this study 
are not highly tapered, all flutter calculations presented herein were 
made by using the uncoupled mode shapes of a uniform cantilever beam. In 
all cases the first torsion mode and the first and second bending modes 
were used. 

Unless otherwise indicated the subsequent discussion deals entirely 
with results obtained with circulation functions calculated by using theo- 
retical two-dimensional lift-curve slopes (eq. ( 7 )) as in references 1 and 

15 . 


Wing designation .- The three-digit system used to identify the wings 
with taper ratio of 0.6 is the same as that used in references 1, 5> and. 
15 . The first digit in this system is the aspect ratio of the full wing 
to the nearest integer. The second and third digits give the quarter- 
chord sweep angle to the nearest degree. For example, wing 445 das an 
aspect ratio of 4, a sweep angle of 45°, and a full -wing taper ratio of 
0.6. Since some of the wings discussed in this paper have identical plan 
forms but different center -of -gravity positions (ref. 6), a single letter 
is appended to the plan-form designation to signify a forward or rearward 
shift in center of gravity. For example, wing 445 das a center of gravity 
at approximately 46 percent chord, whereas the center of gravity of wing 
445F is at about 54 percent chord, and that of wing 445R is at about 
58 percent chord. Wing 400 has a center of gravity at approximately 
45 percent chord, but wing 400R has a center of gravity at about 59 per- 
cent chord. 


\ji a\ -r tr 1 
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For the wing with taper ratio of 1.0, the same system is used, 
except that a fourth digit 1 is added to distinguish the taper ratio . 

Thus wing 4001 has a full-wing aspect ratio of 4, a sweep angle of 0, 
and a taper ratio of 1 . 0 . 

Flutter characteristics .- Calculated flutter characteristics V^V R 
and co^o 3 a are compared in figures 10 to 21 with experimental flutter data 

and with flutter characteristics calculated for the same wings in refer- 
ences 1 and 15 from theoretical distributions of aerodynamic parameters . 
The experimental flutter points shown were obtained at various values of 
density p, whereas, for a particular wing, all of the points calculated 
in this investigation and in reference 1 were obtained at a constant value 
of p which represented approximately an average of the experimental den- 
sities . For each experimental point, however, the normalizing V R was 

calculated by using the appropriate experimental density. On the basis of 
previous experience, it appears that normalizing the experimental flutter 
speeds in this manner accounts for the major effects of density so that 
the resulting ( eX p is considered to depend only slightly on p, at 

least over the range of density variation which occurs herein. The 
variation of V^V R with density is examined in detail in reference 15 . 


The reference flutter speeds V R used in references 5 and 12 for 

wings 400 and 4001 were calculated by employing only two degrees of free- 
dom (first bending and first torsion) . Since three -degree -of -freedom cal- 
culations yield values of V R which are slightly different from the two- 

degree -of -freedom values, the experimental V^V R values for these two 

. , , . , . , . V R (for two degrees of freedom) 

wings have been mutliplied by the ratio — ^7 — 

V R (for three degrees of freedom) 

so that both calculated and experimental flutter-speed ratios as presented 
herein are normalized by Vr for three degrees of freedom. 


Swept Wings — Wings 445, 445F, 445R 

Flutter speeds .- For wings 445, 445F, and 445R values of V/V R 
(figs. 10 to 12) calculated from experimental distributions of C 7 

a,n 

and ac n (for the undeformed wing in steady flow) (figs. 1 and 2 ) are 
in good agreement with experimental V^V R values at subsonic, transonic, 

and supersonic speeds and are also in good agreement at subsonic and 
supersonic speeds with values calculated from theoretical distributions 


14 


of C 7 _ _ and ac n (refs. 1 and 15 )-^. When comparing flutter charac- 


u a,n 


ter is tics obtained from these three sources, however, it should be 


remembered that the experimental distributions of Cj 


a.n 


and ac 


n 


used were obtained from data for a 6-percent-thick wing while the experi- 
mental flutter points were obtained with 4-percent-thick wings. It is 
observed that the dip which appears in the V/V-p curves (figs. 10 to 12) 


in the range 0.8 < M < 1.0 is coincident with a relatively sharp for- 
ward shift in ac n position for the 6-percent-thick wing (fig. 3) • 

Since this ac n shift might be expected to be somewhat less severe for 

a 4-percent-thick wing, use of experimental aerodynamic parameters for a 
4-percent-thick wing might be expected to result in a somewhat shallower 
dip in the V/Vr curves. The above-mentioned effect of ac n change is 


believed to be the principal inaccuracy resulting from use of aerodynamic 

parameters for a 6-percent-thick wing. Figure 3 shows that the values of 

C 7 and ac n for M = 0.9^- and M = O .96 differ significantly from 
*'a,n 

the curves faired through the remainder of the data. Because of this 
deviation the flutter points calculated for these two Mach numbers have 
been ignored in fairing the V/Vp curves in figures 10 to 12. 


Figure 10 shows flutter speeds for wing 445 at supersonic Mach num- 
bers calculated from experimental distributions of C 7 and ac n 

a,n 

obtained both by directly applying equations (l) and (2) to data for 
streamwise sections (solid curves of figs. 1 and 2) and by adapting 
these values near the tip to the strip orientation employed in the flutter 
calculation procedure as described previously (dashed curves of figs. 1 
and 2) . The flutter speeds resulting from these two methods of evalu- 
ating C 7 and ac n are so close together that use of the strip- 

*'a,n 

orientation adaptation is considered unnecessary at least for the wings 
of this investigation. Accordingly, for wings 445F and 445R equations (l) 
and (2) are applied directly for all Mach numbers. 

Figure 10 also shows flutter speeds at subsonic Mach numbers calcu- 
lated by using circulation functions for incompressible flow (C=Fj+iGj). 

It appears that subsonic flutter speeds calculated in this manner for 
wing 445 are low by no more than about 5 percent up to M = 0.90. This 
result further confirms the statement in reference 1 that the modified 


l-For the subsonic points calculated from theoretical aerodynamic 
parameters the value of ac n has been obtained from equation (2) with 

ac = i as discussed in reference 15. The V/Vr and o)/a) a curves shown 
in figures 10 to 21 are identical to those of reference 15. 


^ 4-vo in 


15 


circulation function C 



need be employed only at high 


subsonic and supersonic speeds. Reference 15 discusses this point 

further in connection with flutter calculations employing theoretical 

distributions of C7 and ac n . 

6 a,n n 


Flutter frequencies .- For the three 45° swept wings values of oyo^ 
(figs. 13 to 15) calculated from experimental distributions of C7 

a,n 

and ac n are in good agreement with experimental oyo) a values at sub- 
sonic and transonic speeds and are also in good agreement at subsonic and 
supersonic speeds with values calculated from theoretical distributions 


of C 


and ac n (refs. 1 and 15)* For all three wings at transonic 
speeds the a)/o) a curves calculated from experimental aerodynamic param- 


a,n 


eters turn upward tending to follow the trend of the experimental oyo) a 

values. At Mach numbers past 1.0, however, the upward trend is halted, 
and the curves then tend to follow those calculated from theoretical 
aerodynamic parameters. 


Figure 13 shows flutter frequencies for wing 445 at supersonic Mach 
numbers calculated from experimental distributions of C^ a n and ac n 

obtained both by directly applying equations (l) and (2) to data for 
streamwise .sections (solid curves of figs. 1 and 2) and by adapting 
these values near the tip to the strip orientation employed in the flutter 
calculation procedure (dashed curves of figs. 1 and 2). The closeness of 
the frequencies which result from these two methods of evaluating C^ a n 

and ac n confirms the statement (made previously in connection with 
flutter speeds) that the strip-orientation adaptation is unnecessary here. 


Figure 13 also shows flutter frequencies at subsonic Mach numbers 
calculated by using circulation functions for incompressible flow 
(C = Fj + iGj). The differences between these frequencies and those cal- 


culated by using the modified circulation function 
very small for all subsonic Mach numbers. 


C 



- “1) 


are 


Unswept Wings with Taper Ratio of 0.6 — Wings 400, 400R 

Flutter speeds .- For wings 400 and 400R values of V/V R (figs. l6 

and 17) at M = J~2, calculated from distributions of Ci and acr, 
v 6 a,n n 
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obtained from flight tests of the X-1E airplane, are in much better 
agreement with the experimental flutter data than are the values calcu- 
lated from linearized-theory aerodynamic parameters, especially for 


wing 400. The distributions of 


a,n 


and ac n obtained from the X-1E 


data are compared with values calculated from steady-flow linearized 
theory in figure 6. This comparison shows that the 


<x,n 


distribution 


is accurately predicted by linearized theory except near the fuselage. 
However, as might be expected the linearized theory predicts aerodynamic- 
center positions that are too far rearward. References 1 and 15 hypothe- 
sized that for wing ^00 at supersonic Mach numbers the poor agreement 
between experimental flutter speeds and those calculated from steady-flow 
linearized-theory aerodynamics was related to the close proximity of the 
local aerodynamic centers to the local centers of gravity and to the 
fact that linear theory predicts too-far rearward aerodynamic centers. 

For wing U00 at supersonic Mach numbers the great sensitivity of calcu- 
lated flutter speed to small changes in aerodynamic -center position was 
demonstrated in reference 15 . Comparison of the flutter speeds calcu- 
lated from the X-1E data and from linearized-theory aerodynamic param- 
eters gives further evidence of this sensitivity. In view of this 
sensitivity and in view of the fact that the differences in aerodynamic 
center shown in figure 6 are caused by finite wing thickness, viscosity, 
and other nonlinear effects, it seems doubtful that accurate supersonic 
flutter speeds for this wing could be obtained by using any linearized 
theory . 


For wings 400 and 400R, values of V/V R (figs. l 6 and 17 ) calcu- 
lated from C7 and ac n distributions obtained from the data of 
a,n 

reference 7 (fig. k) appear to be consistent with the experimental 
flutter points for Mach numbers up to 0.8. Since Mj.p « M for these 

unswept wings, and since the circulation function Fq approaches 0 as 
Mjjg approaches 1, the Mach number range near M « M Ti ^ = 1.0 for these 

wings is not accessible to the calculation procedure as given in refer- 
ence 1 . (See previous discussion of this point.) As shown in fig- 
ures l 6 and 17 for Mach numbers near 1.0, the calculated flutter speeds 
appear to increase without limit. 


As indicated previously, an attempt was made to eliminate this 
spurious behavior of the V/Vr curves by employing experimental values 

of the two-dimensional lift-curve slope (fig. 9) in the calculation of 
the circulation function Fq for Mach numbers near 1.0 (Mj^ near 1.0). 
The V/V R values obtained in this manner are shown in figures l6 and 17 . 
It may be seen that although finite VyV-^ values are obtained for 
M = 1.0, the agreement with experimental values still is not good. 


ui o\ ir 1 
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Accordingly, the artifice employed in the computation of the circulation 
functions for these flutter calculations is not used further. 

Values of V/V R for wings 400 and 400R calculated by using total 

lift-curve slope and aerodynamic center (fig. 7) as described previously 
are shown in figures 16 and 17. These values are in good agreement with 
the experimental values up to M = 0.8. Above this Mach number the 
limitation on the circulation function just discussed is again indicated 
by a spurious rise in the calculated flutter speed. 

Flutter frequencies .- For wings 400 and 400R values of co/o^ 

(figs. 18 and 19) at M = v/2~, calculated from distributions of C 7 

t a,n 

and ac n obtained from flight tests of the X-1E airplane (fig. 6), are 
in good agreement with the experimental flutter frequencies. Comparing 
these calculated frequencies with those calculated from linearized- theory 
aerodynamic parameters (refs. 1 and 15) shows that use of aerodynamic 
parameters for the X-1E airplane results in a very large frequency 
reduction for wing 400 but only an insignificant reduction for wing 400R. 
This large frequency change for wing 400 again indicates the great sensi- 
tivity of calculated supersonic flutter characteristics for that wing to 
small changes in aerodynamic-center position. (See previous discussion 
of flutter speeds and reference 15.) Flutter frequencies for wings 400 
and 400R calculated by use of steady-flow aerodynamic parameters obtained 
from the data of reference 7 (f ig. are somewhat closer to the experi- 
mental flutter frequencies than are the curves faired through points 
calculated from linear-theory aerodynamic parameters. Flutter frequencies 
for these two wings calculated by use of experimental values of total 
aerodynamic parameters (fig. 7) closely follow the values calculated from 
distributions of linear-theory aerodynamic parameters . 


Unswept Wing with Taper Ratio of 1.0 — Wing 4001 

Flutter speeds .- Values of V/V R for wing 4001 calculated by using 

total lift-curve slope and aerodynamic center (fig. 8) are shown in fig- 
ure 20. These values are in fairly good agreement with the experimental 
flutter points up to about M = O.85, but above this Mach number a spur- 
ious asymptotic rise in the calculated flutter speed again appears as 
M ~ M tr approaches 1.0. It should be remembered that the aerodynamic 

parameters of figure 8 were obtained with wings with NACA 63AOO4 and 
NACA 65AOO4 airfoil sections (rounded leading edge), while the experi- 
mental flutter points were obtained with wings with a 4-percent-thick 
hexagonal airfoil (sharp leading edge). It is recognized that even 
though the wings of figure 8 and reference 12 have the same thickness, 
these differences of airfoil shape (especially leading-edge radius) can 




18 


lead to significant differences in the lift-curve slope and aerodynamic 
center, particularly at Mach numbers near 1.0. At high subsonic speeds, 
for example, it is probable that the aerodynamic center for the flutter 
wings was slightly farther rearward than the positions indicated in 
figure 8 . If properly accounted for this condition would result in cal- 
culated flutter speeds slightly higher than those shown in figure 20 . 

Flutter frequencies .- Flutter frequencies (fig. 21) calculated by 
using the aerodynamic parameters of figure 8 are slightly lower than 
frequencies calculated by using linear-theory aerodynamic parameters, 
particularly at high subsonic and transonic speeds. 


CONCLUDING REMARKS 


Flutter calculations for several swept and unswept wings through the 
transonic speed range have been made by the modified strip-theory method 
of NACA RM L 57 LIO, which employs steady-flow aerodynamic parameters for 
the undeformed wing. Experimentally determined distributions of steady- 
flow aerodynamic parameters for the undeformed wings were used in the 
present calculations . Comparisons of these calculated results with 
experimental flutter data and with calculations previously made by using 
linearized-theory static aerodynamic parameters indicate that the method 
employed gives accurate flutter results for swept wings at subsonic, 
transonic, and supersonic speeds. However, since this method of flutter 
calculation is not applicable when the Mach number component normal to 
the leading edge is near 1 . 0 , it appears that the transonic flutter 
characteristics of unswept wings cannot be calculated by the method as 
given in NACA RM L 57 LIO. An attempt was made to remove this limitation 
by using experimental values of two-dimensional lift-curve slope (instead 
of theoretical values) in the calculation of circulation functions when 
the Mach number component normal to the leading edge was near 1.0. 

Applying this procedure to two unswept wings removed the spurious asymp- 
totic rise of flutter speed near Mach number 1.0, but it did not result 
in very close agreement between calculated and experimental flutter speeds . 

The possibility of obtaining estimates of flutter characteristics for 
unswept wings by use of total aerodynamic parameters has also been exam- 
ined. Although no supersonic calculations were made using total aerodyna- 
mic parameters, it appears that the procedure outlined gives reasonable 
estimates of flutter characteristics at subsonic speeds. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., August 26, 1959* 
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TABLE I.- SUMMARY OF CONFIGURATIONS FOR WHICH FLUTTER CALCULATIONS WERE MADE INCLUDING 


SOURCES OF THE EXPERIMENTAL STATIC AERODYNAMIC PARAMETERS USED 



Configuration for which flutter tests 
and flutter calculations were made 

Source of experimental static 
aerodynamic parameters 

Use of 

experimental 

aerodynamic 

parameters 

Flutter 
results 
presented 
herein 
in figure 

Wing 

designation 

Aspect 

ratio 

Sweep 
angle 
(quarter- 
chord) , 
deg 

Taper 

ratio 

Airfoil 

section 

(stream- 

wise) 

Center of 
gravity, 
percent 
chord 

Reference 
in which 
presented 

Figure 

herein 

Aspect 

ratio 

Sweep 

angle, 

deg 

Taper 

ratio 

Airfoil 

section 

(stream- 

wise) 

44^ 

4.0 

45 

0.6 

NACA 

65A004 

46 

3, 4 

1 , 2 , 3 

4.0 

A cA = 1 * 5 

0.6 

NACA 

65 AOO 6 

Spanwise 

distributions 

. 10, 13 

445F 





34 

3 , 4 

l, 2 , 3 

4.0 

a cA = 45 

.6 

— do— 

do 

11 , 14 

445R 





58 

3, 4 

1 , 2 , 3 

4.0 

A cA = 45 

.6 

—do— 

do 

12 , 15 

400 

4.0 

0 

0.6 

NACA 

65A004 

45 

7 

4, 5 

4.0 

a c /2 = 0 

0.5 

NACA 

65A004 

Spanwise 

distributions 

16 , 18 

400R 





59 

7 

4, 5 

4.0 

A c /2 = 0 

• 5 

— do— 

do 

17 , 19 

400 





45 

Unpublished 

6 

4.0 

A .4c “ 0 

•5 

NACA 

64A004 

do 

16 , 18 

400R 





59 

Unpublished 

6 

4.0 

A .4c = 0 

•5 

— do — 

do 

17 , 19 

400 





45 

9 

7(a) 

4.0 

\/k “ 0 

.6 

NACA 

65A004 

Total values 

16 , 18 

400R 





59 

9 

7(a) 

4.0 

A c/k = 0 

.6 

— do— 

do 

17 , 19 

400 





45 

14 

9 

--- 

— 

--- 

—do— 

Two-dimensional 

C 2 in F + iG 
c a 

16 , 18 

400R 





59 

14 

9 

— 

— 

— 

— do — 

do 

17 , 19 

4001 

4.0 

0 

1.0 

Hexagon 

50 

10 

8(a) 

4.0 

A = 0 

1.0 

NACA 

65A004 

Total values 

20, 21 






50 

11 

8(a) 

4.0 

A = 0 

1.0 

NACA 

63A004 

do 

20, 21 


ro 
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(a) M = 0.60. 



0 .1 .2 


• 3 *U .5 -6 .7 

*1 


.8 .9 1.0 


(b) M = 0.80. 

Figure 1.- Distributions of steady-flow aerodynamic parameters obtained 
from the measured load distributions of reference 3 for a wing of 
A = 4.0, A = 0.6, A c y^ = 45°, and NACA 65AOO6 airfoil. Symbols 
indicate values used in flutter calculations. Solid curves were 
obtained directly from equations (l) and (2). Dash curves were 
obtained by altering the results from equations (l) and (2) to con- 
form to strip theory near the tip. 
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(c) M = 0.90. 



(d) M = 0.93. 


Figure 1.- Continued. 
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(e) M = 0.96. 



r\ 


(f) M = 1.20. 


Figure 1.- Concluded. 
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(a) M = 0.94. 



Ob) M = 0.97. 

Figure 2.- Distributions of steady-flow aerodynamic parameters obtained 
from the measured load distributions of reference 4 for a wing of 
A * 4.0, A = 0.6, Ac = 45°, and NACA 65AOO6 airfoil. Symbols 
indicate values used in flutter calculations . Solid curves were 
obtained directly from equations (l) and (2). Dash curves were 
obtained by altering the results from equations (l) and (2) to con- 
form to strip theory near the tip. 
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(c) M = 0.99. 



(d) M = 1.02. 


Figure 2.- Continued. 
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(e) M = 1.11. 



(f) M = 1.13. 


Figure 2.- Concluded. 



Figure 3 • - Variation of steady-flow aerodynamic parameters with Mach number for several spanwise 
stations on the wing of reference 3 (plain symbols) and reference 4 (flagged symbols). 

A = 4.0; A = 0.6; Ac /4 = 45°; NACA 65 AOO 6 airfoil. 
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o .1 .2 .3 Ji .5 .6 .? .8 .9 l.o 
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(a) M = 0.60. 





(b) M = 0.70. 

Figure 4.- Distribution of steady-flow aerodynamic parameters obtained 
from the measured pressure distributions of reference 7 for a wing 
of A = 4.0, A = 0.5, A c / 2 = 0, and NACA 65A004 airfoil. Symbols 
indicate values used in flutter calculations. 
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(c) M = 0.80. 


c la,n 



•1 .2 .3 .1* .5 .6 .7 .8 .9 1.0 


(d) M = 0.85. 
Figure 4.- Continued. 


31 



-.8 


-.7 


-.6 


-.5 


ac n 


:h 


-.5 


-.2 


-.1 


0 


0 .1 .2 .3 .1+ .5 .6 .7 .8 .9 l.o 




(e) M = 0.88. 



(f) M = 0.90. 
Figure 4.- Continued. 
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(g) M = 0.92. 



I 


(h) M = 0.94. 
Figure 4.- Continued. 
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(i) M = 0.96. 



(j) M = 0.98. 
Figure 4.- Continued. 
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(k) M = 1.00. 



( l ) M = 1.05. 
Figure 4.- Concluded. 
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Figure 5*~ Variation of steady-flow aerodynamic parameters with Mach 
number for three spanwise stations on the wing of reference 7 . 

A = 4.0; A = 0.5; Ac /2 = 0; NACA 65A004 airfoil. 
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-b 


ac n 


• 5 


-.2 


-.1 
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T) 

Figure 6.- Distributions of steady-flow aerodynamic parameters obtained from flight test data 
for the X-1E airplane at M = 1.4l. A = 4.0; A = 0.5; = 0; modified NACA 64A004 air- 

foil. Symbols indicate values used in flutter calculations. 
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(a) Experimental steady-flow total aerodynamic parameters obtained 
from reference 9 for a wing of A = 4.0, A = 0.6, = 0, 

and NACA 6^A00k airfoil. ' 



•1 *2 .3 .1* .5 .6 .7 .8 .9 l.o 


(b) Calculated spanwise distribution of Cz a ^ n which was used in 

combination with the experimental steady-flow total aerodynamic 
parameters of figure 7 (&) in flutter calculations for wings 400 
and 400R. Values of ^Z a ^ n shown were calculated by the method 

of reference 8 for M = O. 75 . 

Figure 7*~ Aerodynamic quantities involved in the estimation of 
flutter characteristics for wings kOO and kOOR by use of total 
aerodynamic parameters. 
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(a) Experimental steady-flow total aerodynamic parameters for wings 
of A 4.0, A = 1.0, A = 0. Plain symbols represent data from 
reference 11 (wing with NACA 65A004 airfoil), and flagged symbols 
represent data from reference 10 (wing with NACA 65A004 airfoil). 



0 .1 .2 .3 .U .5 .6 .7 .8 .9 1.0 

■n 


(b) Calculated spanwise distribution of C? which was used in 

QLy n 

combination with the experimental steady-flow total aerodynamic 
parameters of figure 8 (a) in flutter calculations for wing 4001. 
Values of Cj a ^ n shown were calculated by the method of refer- 
ence 8 for M = O. 75 . 

Figure 8 .- Aerodynamic quantities involved in the estimation of 
flutter characteristics for wing 4001 by use of total aero- 
dynamic parameters. 
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Figure 9 .- Experimental two-dimensional steady-flow lift-curve slope 
for NACA 65A004 airfoil at zero lift as given in reference 14 and 
the relation of this quantity to linear-theory two-dimensional 
steady-flow lift-curve slope. 
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Figure 10.- Variation of flutter speed with Mach number for wing For calculated points 

p = 0.003800 slug/cu ft and Vp = 735-0 ft /sec. 


O Flutter experiment 
Flutter calculation using 

O Theoretical Cj a and ac n from reference 1 

A Experimental Cj a n ai d ac n from figures 1 and 2 

d Experimental Cj a n and ac n altered to conform to 
strip theory near the tip as on figures 1 and 2 
Experimental ^ and ac n from figures 1 and 2 
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Figure 11.- Variation of flutter speed with Mach number for wing 445F. For calculated points 

p = 0.003000 slug/cu ft and V R = 840.0 ft /sec. 
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Figure 12.- Variation of flutter speed with Mach number for wing 44 5R. For calculated points 

p = O.OO2378 slug/cu ft and Vp = 928.7 ft/sec. 
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O Flutter experiment 

Flutter calculation using 

□ Theoretical Cj a ^ n and ac n from reference i 
A Experimental Cj a>n and ac n from figures 1 and 2 
^ Experimental Cj a ^ n and ac n altered to conform to 

strip theory near the tip as on figures 1 and 2 
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Figure 13.- Variation of flutter frequency with Mach number for wing 44-5. For calculated points 

p = 0.003800 slug/cu ft and = 2,192 radians/sec. 
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O Flutter experiment 
Flutter calculation using 

□ Theoretical Cj a n and ac n from reference 1 

A Experimental Cj and ac n from figures 1 and 2 




















0 .1 .2 .3 . 1 + .5 -6 .7 .8 .9 1.0 l.l 1.2 1.3 l . i + 

M 


P' 

pr 


Figure 14.- Variation of flutter frequency with Mach number for wing 445F. For calculated points 

p = 0.003000 slug/cu ft and = 1,144 radians/sec. 
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Figure 15.- Variation of flutter frequency with Mach number for wing 445 R. For calculated points 

p = O.OO2378 slug/cu ft and = 2,306 radians/sec. 
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Figure l6.- Variation of flutter speed with Mach number for 
wing 400. For calculated points p = 0.002378 slug/cu ft 
and Vr = 976.5 ft/sec. 
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Figure 17- - Variation of flutter speed with Mach number for wing kOOR. For calculated points 

p = 0.003100 slug/cu ft and V R = 852.5 ft /sec. 
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Figure 18.- Variation of flutter frequency with Mach number for wing 4-00. For calculated points 

p = 0.002378 slug/cu ft and cD a = 2,463 radians/sec . 
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Figure 19*- Variation of flutter frequency with Mach number for wing 4-00R. For calculated points 

p = 0.003100 slug/cu ft and = 1,982 radians/sec. 
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Figure 20.- Variation of flutter speed with Mach number for 
wing 4001. For calculated points p = 0.002378 slug/cu ft 
and Vr = 828.5 ft/sec. 
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Figure 21.- Variation of flutter frequency with Mach number for wing 4001. For calculated points 

p = 0.002578 slug/cu ft and o^, = 2,048 radians/sec. 


VJ 1 

H 





